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Key Points 


Box model simulated glyoxal production from three isoprene oxidation 
mechanisms differ greatly 

Aerosol uptake of glyoxal was constrained using airborne in situ measurements 
and a global model 

Model results show that glyoxal contributes 0-14% of SOA in the Southeast U.S. 
during summer 
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Abstract 


We use a 0-D photochemical box model and a 3-D global chemistry-climate model, 
combined with observations from the NOAA Southeast Nexus (SENEX) aircraft 
campaign, to understand the sources and sinks of glyoxal over the Southeast United 
States. Box model simulations suggest a large difference in glyoxal production among 
three isoprene oxidation mechanisms (AM3ST, AM3B, and MCM v3.3.1). These 
mechanisms are then implemented into a 3-D global chemistry-climate model. 
Comparison with field observations shows that the average vertical profile of glyoxal is 
best reproduced by AM3ST with an effective reactive uptake coefficient Ye1yx of 2 x 10°, 
and AM3B without heterogeneous loss of glyoxal. The two mechanisms lead to 0-0.8 ug 
m° secondary organic aerosol (SOA) from glyoxal in the boundary layer of the Southeast 
U.S. in summer. We consider this to be the lower limit for the contribution of glyoxal to 
SOA, as other sources of glyoxal other than isoprene are not included in our model. In 
addition, we find that AM3B shows better agreement on both formaldehyde and the 
correlation between glyoxal and formaldehyde (Rgg= [GLYX]/[HCHO)]), resulting from 
the suppression of 5-isoprene peroxy radicals (6-ISOPO2). We also find that MCM v3.3.1 
may underestimate glyoxal production from isoprene oxidation, in part due to an 
underestimated yield from the reaction of IEPOX peroxy radicals (IEPOXOO) with HO>. 
Our work highlights that the gas-phase production of glyoxal represents a large 


uncertainty in quantifying its contribution to SOA. 


1 Introduction 


Glyoxal (CHOCHO) is one of the most abundant dicarbonyl compounds in the 


atmosphere. Its sources include direct emissions from biofuel use and biomass burning, 
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and secondary production from oxidation of various volatile organic compounds (VOCs) 
[Fu et al., 2008; Hays et al., 2002; Myriokefalitakis et al., 2008]. Glyoxal has a lifetime 
of about 1-3 hours against photolysis and oxidation by OH at midday [Feierabend et al., 
2009; Volkamer et al., 2005; Washenfelder et al., 2011]. It is highly water-soluble, with a 
Henry’s Law constant of 3.0 - 4.2 x 10° M atm! at 298 K [Sander, 2015]. In aerosol 
water, the solubility increases rapidly at low salt concentrations (up to 5.0 x 10° M atm’) 
[Kampf et al., 2013; Waxman et al., 2015], due to the formation of glyoxal hydrate - 
sulfate complexes (“salting-in’”) [Kurtén et al., 2015]. However, this increase is inhibited 
at high salt concentrations by the kinetic limitation of gas-particle partitioning [Kampf et 
al., 2013; Knote et al., 2014]. Laboratory and field studies showed that glyoxal readily 
undergoes heterogeneous uptake to aerosols and cloud droplets to form secondary 
organic aerosol (SOA) [Carlton et al., 2007; Ervens et al., 2011; Galloway et al., 2009; 


Galloway et al., 2011b; Liggio et al., 2005a; b; Lim et al., 2005; Volkamer et al., 2007]. 


Glyoxal can provide important constraints on quantifying VOC emissions and 
oxidation mechanisms. Observations from aircraft, ground and satellite show that glyoxal 
is often highly correlated with formaldehyde (HCHO), another product of VOC oxidation 
[DiGangi et al., 2012; Kaiser et al., 2015; MacDonald et al., 2012; Miller et al., 2014; 
Stavrakou et al., 2009; Vrekoussis et al., 2010]. The ratio of glyoxal to formaldehyde 
(surface or tropospheric column concentrations), Rg, varies for different biogenic and 
anthropogenic VOC precursors, so it can be used to quantify the source strength of local 
VOC emissions [DiGangi et al., 2012; Kaiser et al., 2015; Miller et al., 2014; Vrekoussis 
et al., 2010]. The NOx (NO, = NO + NO») dependence of Rep, although it varies greatly 


for different VOCs and may offer additional information on VOC oxidation [DiGangi et 
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al., 2012; Kaiser et al., 2015; Vrekoussis et al., 2010]. Thus, it is important to evaluate 
model performance on Rgr and its NOx dependence on regional and global scales, 
providing critical information for present and future satellite validation [Vrekoussis et al., 


2010]. 


Sources and sinks of glyoxal remain largely uncertain [Ervens et al., 2011; Fu et al., 
2008; Galloway et al., 2011a; Myriokefalitakis et al., 2008; Stavrakou et al., 2009; 
Volkamer et al., 2007; Washenfelder et al., 2011]. Production from isoprene oxidation 
represents a major source of glyoxal on the global scale [Fu et al., 2008]. Chamber 
experiments suggest a significant yield (up to 3%) of glyoxal from the first generation of 
isoprene oxidation under high NO, conditions (500 ppbv NO) [Galloway et al., 201 1a; 
Volkamer et al., 2006]. As we show below, this high yield may not reflect the distribution 
of B- and 6-isoprene peroxy radicals (ISOPOz2) as a function of their lifetimes under 
ambient conditions [Peeters et al., 2014]. Another large uncertainty lies in the 
heterogeneous loss of glyoxal to aerosols and cloud droplets, which contributes to SOA 
mass. Laboratory studies indicate that the uptake of glyoxal depends on aerosol and cloud 
composition [Corrigan et al., 2008; Jang et al., 2002; Kroll et al., 2005; Liggio et al., 
2005a; b; Noziére et al., 2009; Volkamer et al., 2009; Waxman et al., 2015], ambient 
relative humidity (RH) [Corrigan et al., 2008; Hastings et al., 2005; Kampf et al., 2013; 
Liggio et al., 2005a] and temperature [Gomez et al., 2015]. By assuming irreversible 
reactive uptake with an uptake coefficient Ye1yx = 2.9 x 10°, Fu et al. [2008] found that 
this process accounts for 14% of glyoxal loss globally. Volkamer et al. [2007] estimated 
Yelyx Of 3.7 Xx 10° in Mexico City, where glyoxal contributes about 15% of ambient SOA. 


A much lower Yelyx (0-8 x 10“ for the day and (2 + 1) x 10“ for the night) was derived in 
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Los Angeles [Washenfelder et al., 2011], where kinetic limitations at the high salt 
concentrations likely reduce the rate of SOA formation [Kampf et al., 2013; Knote et al., 


2014]. 


The Southeast Nexus (SENEX) aircraft campaign, which took place in June—July of 
2013, aimed to improve the understanding of the interactions between biogenic and 
anthropogenic emissions over the Southeast U.S. It provided a detailed characterization 
of tropospheric photochemistry (gas and aerosol), including ozone, NO, (NO, NO> and its 
atmospheric oxidation products), biogenic VOCs (isoprene, terpene), isoprene oxidation 
products and organic aerosols [Warneke et al., 2016]. In particular, glyoxal was measured 
on board the NOAA WP-3D aircraft using a cavity enhanced absorption spectroscopy 
technique [Min et al., 2015], and formaldehyde was measured using laser induced 
fluorescence technique [Cazorla et al., 2015]. These measurements provide an 
unprecedented opportunity to examine the sources and sinks of glyoxal in this region, as 


well as its contribution to SOA mass. 


Here we first examine the glyoxal yield from isoprene oxidation in a box model with 
three chemical mechanisms (AM3ST, AM3B, and the Master Chemical Mechanism 
v3.3.1 [Jenkin et al., 2015]). We then use field observations from the SENEX field 
campaign, interpreted with a chemistry-climate model, to understand the sources and 
sinks of glyoxal over the Southeast U.S. The comparison of the model to high resolution 
aircraft data for both glyoxal and formaldehyde provides important new constraints on 
the potential for glyoxal to form SOA, but also highlights uncertainty in the mechanism 


for isoprene oxidation, the single largest source of glyoxal in the Southeast U.S. 
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2 Methods 


AM3 is the atmospheric component of the Geophysical Fluid Dynamics 
Laboratory (GFDL) coupled model CM3. The dynamical core, physical 
parameterizations, cloud and precipitation processes, and cloud-aerosol interactions in 
AMs3 are described in detail in Donner et al. [2011]. Chemistry in a previous version of 
AM3 has been described by Naik et al. [2013]. In this work, we nudge the horizontal 
wind field in the model toward values from the NCEP Global Forecast System (GES), 
allowing the model to simulate synoptic conditions corresponding to those sampled 
during field campaigns [Lin et al., 2012]. We also apply finer vertical grids for 
convection plumes than the standard AM3 to improve the wet removal of tracers during 
summer time [Paulot et al., 2015]. The photolysis module has been updated to Fast-JX 
v7.1 (ftp://128.200.14.8/public/prather/Fast-J), to compute the impacts from aerosols and 
clouds interactively. Dry deposition velocities prescribed in the model reflect rapid dry 


deposition of oxidized organic compounds [Nguyen et al., 2015]. 


Biogenic emission of isoprene is computed using the Model of Emissions of Gases 
and Aerosols from Nature (MEGAN) inventory [Guenther et al., 2006], with a total 15.9 
Tg C in North America in June-August of 2013, which is slightly higher than previous 
estimates of 12.2-14.6 Tg C (June-August of 2006) in the same region [Millet et al., 
2008]. We reduce the isoprene emissions estimated by MEGAN by 20% (to 12.7 Tg C), 
to be consistent with other estimates of isoprene emission over the Southeast U.S. 
[Warneke et al., 2010]. We do not consider glyoxal production from oxidation of terpenes 
and aromatics in this work, as their contribution is relatively small compared to isoprene 


over the Southeast U.S. [Kaiser et al., 2015]. Inclusion of larger isoprene emissions and 
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other VOC sources would proportionally increase the magnitude of glyoxal sinks derived 
in the analysis below. Anthropogenic emissions in 2013 follow the Representative 
Concentration Pathway 8.5 (RCP 8.5) projection [Lamarque et al., 2011]. We applied 
diurnal variation anthropogenic NO, emissions in North America following Mao et al. 
[2013b]. Anthropogenic NOx emissions in RCP 8.5 are 0.34 Tg N/month over North 
America, comparable to the estimate from NEI 2011 of 0.29 Tg N/month [Travis et al., 
2015]. We reduce these anthropogenic NO, emissions by 25% to 0.26 Tg N/month, to be 


consistent with recent estimates in this region [Anderson et al., 2014]. 


2.1 Isoprene chemistry 


The standard isoprene mechanism used in AM3, “AM3ST”, is largely based on Mao 
et al. [2013b] following chamber observations from Paulot et al. [2009a; 2009b] and 
Crounse et al. [2011]. The main features of this mechanism include: (1) substantial yield 
(7%) of first generation organic nitrates from 6-ISOPO2, compared to 4.7% from B- 
ISOPO2; (2) detailed assignment of the fate of first- and second-generation organic 
nitrates; and (3) isomerization of ISOPO2 using the laboratory-determined rate constants. 
We further update this chemistry in several aspects. First, we assume a 25% molar yield 
of glyoxal from ISOPO 2 isomerization as the first-generation product [Marais et al., 
2016]. Second, the reaction of ISOPO2 with HOz2 is updated following St. Clair et al. 
[2016a] to reflect a higher yield of unsaturated hydroxy hydroperoxides (ISOPOOH) 
(94%). Third, we adopt a new methylvinyl ketone (MVK) oxidation chemistry in which 
the glycolaldehyde molar yield is increased from 0.53 to 0.72 under high NO, conditions 
[Praske et al., 2015]. Fourth, we include fast photolysis of carbonyl organic nitrates 


[Miiller et al., 2014]. Fifth, we adopt a substantially slower ozonolysis rate of isoprene B- 
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hydroxy nitrate (ISOPNB, 3.7 x 10°? molec’'cm’s" [Lee et al., 2014]) than the previous 
value from Lockwood et al. [2010] (1.06 x 10°! molec cms”). Lastly, we update the 
NO3-initiated chemistry of isoprene following MCM v3.2 [Jenkin et al., 1997; Saunders 
et al., 2003]. As illustrated in Figure S1, glyoxal is produced in the first generation of 
isoprene oxidation, both from the decomposition of alkoxyl radical (DIBOO shown in 
Figure $1) through the 6-channel under high NO, conditions [Peeters and Nguyen, 2012], 
and from isomerization of ISOPO,. under low NO, conditions [Stavrakou et al., 2010]. 
Glyoxal is further produced in later generations from oxidation of carbonyl compounds 


including glycolaldehyde and isoprene epoxydiol (IEPOX). 


In this work, we introduce an additional mechanism, “AM3_beta” (AM3B), to reflect 
more recent updates to the understanding of isoprene oxidation. Both theoretical and 
experimental studies indicate that B- and 6-[SOPO, undergo fast interconversion [Bates et 
al., 2014; Crounse et al., 2011; Peeters et al., 2014]. One important implication is that 
oxidation products observed in experimental chambers at ~500 ppbv NO may not reflect 
the products in ambient air where NO is several orders of magnitude lower. This is due to 
a higher fraction of 5-ISOPO, loss via bimolecular reactions than interconversion back to 
B-ISOPO, under high NO, conditions. To reflect this change, we only allow B-ISOPO; to 
react with NO because 5-ISOPO) is believed to only account for < 3% of total ISOPO, 
under ambient conditions [Peeters et al., 2014]. The yield of B-hydroxyl isoprene nitrate 
(ISOPNB) is assumed to be 0.1, an average of the suggested values 0.09 [Xiong et al., 
2015] and 0.16 [Teng et al., 2015]. The yields of HCHO, MVK and methacrolein 


(MACR) from B-ISOPO, + NO are adjusted accordingly for carbon balance. As we show 
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below, this update changes the production of formaldehyde and glyoxal under high NOx 


conditions. 


The Leeds Master Chemical Mechanism (MCM) is a near-explicit chemical 
mechanism that describes gas-phase VOC chemistry. The latest revision, v3.3.1, includes 
recent updates for OH-initiated isoprene chemistry [Jenkin et al., 2015], such as isoprene- 
peroxy radical interconversion and isomerization pathways [Peeters et al., 2014], and the 
chemistry of IEPOX [Bates et al., 2014; Paulot et al., 2009b]. Consequently, these 
updates show significant differences in simulated HO, (HO, = OH + HO 2), NO, and 
major oxidation products compared to earlier versions of MCM. To evaluate the 
performance of MCM v3.3.1 in our global model, we implement a MCM-like mechanism 
(‘AM3M_’) by adjusting the production of glyoxal from major pathways in the AM3B 
mechanism, to approximate the glyoxal and HCHO yields of MCM v3.3.1 in a highly 


condensed chemical mechanism suitable for use in a global model (Figure S2). 
2.2 Heterogeneous loss of glyoxal and methylglyoxal 


Heterogeneous loss onto aerosols and cloud droplets plays an important role in the 
fate of dicarbonyls. Here we assume this process is irreversible and can be represented by 


a first-order reactive uptake rate constant k [Mao et al., 2010] 
ba a A 
D, vy 
where @ is the effective radius of aerosols or cloud droplets, D, is the gas phase 


diffusion constant, 7 is the reactive uptake coefficient, A is the surface area of aerosols 


or cloud droplets, and v is the mean molecular velocity of the gas molecule. 
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We allow this process to take place on five types of aerosols including sulfate, black 
carbon (BC), primary organic carbon (POC), sea salt, mineral dust and SOA, with 
hygroscopic growth included [Mao et al., 2013a]. Given the large uncertainties associated 
with the uptake of glyoxal, we use field observations to assess the possible range of Ygiyx. 
We do not include glyoxal uptake by cloud droplets in this work, because the uptake 
coefficient by aerosols is much greater, and because large uncertainties are associated 
with mixing of cloudy and no-cloud volumes within a grid box [Huijnen et al., 2014; 
Jacob, 2000]. We include heterogeneous loss of methylglyoxal with an uptake coefficient 
of 7.6 x 10° [Zhao et al., 2006], to produce better closure with total organic aerosol (OA) 


in the model. 


The irreversible surface uptake coefficient, Yglyx applied here should be viewed as an 
“effective reactive uptake coefficient”, to simply represent the net heterogeneous loss of 
glyoxal that does not revert back to the gas-phase, and provide an estimate for its 
contribution to SOA mass on the regional scale. Previous studies have indicated that this 
process is likely reversible to some extent [Galloway et al., 2009; Kroll et al., 2005]. We 


show below that the values of Y¢1yx are dependent on the choice of gas-phase chemistry. 
3 NO,-Dependent Glyoxal and HCHO Yields 


We first conduct box model simulations to quantify the yield of glyoxal and HCHO 
from isoprene oxidation as a function of NO, levels, using the Dynamically Simple 
Model of Atmospheric Chemical Complexity (DSMACC) [Emmerson and Evans, 2009]. 
We test three mechanisms (AM3ST, AM3B, and MCM v3.3.1) in the model. The model 


is initialized at 8:00 local time (LT) for mid-latitude summer conditions, with 1 ppbv of 
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isoprene, 60 ppbv of ozone, 150 ppbv of CO, and a choice of 0.01, 0.1 or 1 ppbv of NOx. 
While NO, is held constant, isoprene is allowed to decay over time. Figure | shows the 
cumulative yields of glyoxal and HCHO per unit carbon from isoprene oxidation at 
various NO, levels after three days of model integration with their loss processes turned 
off [Palmer et al., 2003]. We also use tagged tracers for individual pathways (f/6- 
TSOPO2 + NO, ISOPO2 + HO, isomerization of ISOPO2) to compute their contribution in 


each mechanism. 


A striking difference between the MCM v3.3.1 and AM3 mechanisms is the glyoxal 
production under low NO, conditions. While the AM3 mechanisms show a large yield of 
glyoxal from isomerization of ISOPO, and the oxidation of IEPOX (via the ISOPO, + 
HO, pathway), MCM v3.3.1 shows very little production of glyoxal from these channels. 
For the isomerization pathway, MCM v3.3.1 assumes 50% of 1,6-H shift isomerization 
flux produces hydroperoxy-aldehydes (HPALDs), with another 50% produces complex 
dihydroperoxy formyl peroxy radicals, referred as “di-HPCARPs” by Peeters et al. 
[2014]. MCM v3.3.1 includes low production of glyoxal from the degradation of 
HPALDs and di-HPCARPs. On the other hand, the AM3 mechanisms do not include di- 
HPCARPs, and assume fast photolysis as the only loss process of HPALDs, with prompt 


yield of glyoxal following Stavrakou et al. [2010]. 


Another large difference in glyoxal yield comes from the treatment of the ISOPO2 + 
HO, pathway, particularly regarding the fate of IEPOX. MCM v3.3.1 assumes a major 
production of organic peroxides from the reaction of IEPOX peroxy radicals (IEPOXOO) 
(a-carbonyl peroxy radicals) with HO, while AM3 mechanisms follow Paulot et al. 


[2009b] and assume full radical propagation for this reaction, with a 28% yield of glyoxal. 
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As a result, a much higher yield of glyoxal is predicted from the ISOPO2 + HO> pathway 
in the AM3 mechanisms compared to MCM v3.3.1. In fact, several studies have 
confirmed significant OH production from the reaction of HO with a-carbonyl peroxy 
radicals [Dillon and Crowley, 2008; Hasson et al., 2004; Jenkin et al., 2007]. This is also 
consistent with recent studies of IEPOX kinetics, which show very little formation of 


peroxides under low NO, conditions [Bates et al., 2014; Jacobs et al., 2013]. 


Under high NO, conditions, the three mechanisms show better agreement for glyoxal 
production, though with different pathways contributing. MCM v3.3.1 shows a dominant 
production from the B-ISOPO, + NO pathway, mainly due to a higher yield of glyoxal 
from oxidation of glycolaldehyde (20%), based on Niki et al. [1987]. An even higher 
yield of glyoxal (29%) was adopted by Galloway et al. [2011a]. In contrast, the AM3 
mechanisms apply a lower yield of glyoxal (13%) following Butkovskaya et al. [2006]. 
Surprisingly, both AM3ST and AM3B show significant production from isomerization of 
ISOPO> and the ISOPO2 + HO» pathway, suggesting the potential importance of these 
channels for glyoxal production even under high NOx conditions. AM3ST has a major 
production from the 5-ISOPO2 + NO pathway, leading to 60% higher glyoxal yield than 
AM3B. MCM v3.3.1 shows a small contribution from the 6-[SOPO2 + NO channel, due 
to a small fraction of 6-[SOPO> in the radical pool of ISOPO2 (< 3%). This reflects a 
different distribution of 5-ISOPO, and B-ISOPOz2 under ambient conditions vs. very high 
NO chamber conditions [Peeters et al., 2014], consistent with our assumption in AM3B. 
Overall, the AM3 mechanisms show decreasing glyoxal production with increasing NO, 


concentration, whereas MCM v3.3.1 shows the opposite. 
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The HCHO yield and its NO, dependence appear to be more consistent across the 
different mechanisms. Under low NO, conditions, isomerization of ISOPO. and the 
ISOPO,2 + HO: pathway contribute most to HCHO. At NO,= 1 ppbv, the ISOPO2 + NO 
pathway becomes the primary source of HCHO. The HCHO yield in AM3B is 17% 
higher than AM3ST and comparable to MCM v3.3.1, mainly due to a higher HCHO yield 
from the B-ISOPO, + NO pathway than the 6-ISOPO2 + NO pathway. It should be noted 
that the cumulative yields shown in Figure 1 represent the potential glyoxal and HCHO 
produced from sufficient oxidation of 1 ppbv isoprene and the intermediate VOCs. This 
is different from ambient conditions that isoprene is continuously emitted and the 
production of glyoxal and HCHO are dependent on OH levels (i.e. the production rate of 
ISOPO2). Chemical loss of HCHO and glyoxal are similar across the mechanisms (Table 


S1). 


The NO, dependent yields of glyoxal and HCHO from box model results are useful to 
identify their major production pathways at different NOx levels. This knowledge can be 
translated to 3D model outputs and help evaluate model performance from these 


pathways using aircraft observations. 


4 Observational Constrains on Glyoxal Production 


4.1 Vertical profiles of gaseous and particulate species 


We compare the AM3 model predictions for glyoxal, formaldehyde, and isoprene 
oxidation products to measurements of these species acquired during the SENEX field 
campaign. All measurements are averaged to | minute time resolution. The measurement 


accuracies for isoprene, formaldehyde, NO,, glyoxal, organic aerosol and aerosol surface 
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area are 25%, 10%, 5%, 5.8%, 50%, and 36% respectively ([Warneke et al., 2016] and 
the references therein). Although other carbonyl compounds have recently been found to 
suffer from inlet artifacts arising from catalytic conversion of organic hydroperoxides on 
metal inlet surfaces [Rivera-Rios et al., 2014], we do not expect such interferences for 
glyoxal since no metal surfaces are present in the glyoxal sampling line [Min et al., 2015], 
and since there are currently no known mechanisms for conversion of hydroperoxides to 
glyoxal. We also expect little interferences for formaldehyde since its residence time 
during exposure to metal surfaces is very small [Cazorla et al., 2015], This is confirmed 
by recent laboratory tests, suggesting the interference is < 5% for this specific instrument 
[St. Clair et al., 2016b]. We exclude biomass burning, urban plumes, stratospheric air 
(CH3CN = 225 pptv, NO,/NO, > 0.4 or NO» > 4 ppbv or O3/CO > 1.25 mol/mol, 
respectively) from our analysis following Hudman et al. [2007], and omit data from the 
Ozark Mountains, where the model shows a significant positive bias for isoprene (Figure 


S3). We also exclude nighttime flights from our analysis. 


Figure 2 shows the mean vertical profiles of observed and modeled isoprene, glyoxal, 
HCHO and other related species during SENEX. Model output is sampled along the 
flight tracks and at the flight time with 1 hour time resolution. We include two model 
simulations (AM3ST and AM3B). MCM v3.3.1 is not included because of its complexity. 
Instead, we show AM3M in our global model as a proxy to MCM v3.3.1 in the following 
sections. Simulated NO, agrees with the observations at all altitudes except for a slight 
overestimate (~20%) near the surface. Observed isoprene concentration peaks near the 
surface at 1.2 ppbv and decreases gradually with altitude. Both AM3ST and AM3B show 


a positive bias of isoprene below | km, but no such bias is evident for HCHO. As HCHO 
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in this region is dominantly produced from isoprene and has a relatively longer lifetime 
than that of isoprene (4 h vs. < 1 h), we attribute this positive bias of modeled isoprene 
partly to sampling bias and partly to shallow boundary layer and slow vertical mixing in 
the current model. We further examine modeled OH and J values using ICARTT (2004) 
[Mao et al., 2013b] and NOMADSS (2013) 
(https://www.eol.ucar.edu/field_projects/nomadss) field observations over the eastern 
U.S., and we find that the model agrees with observed OH and J values within 20% 


(Figure S4 and S5). 


Like isoprene, both HCHO and glyoxal decrease with altitude, reflecting their short 
photochemical lifetimes (4 h and 3.5 h) and their dominant source from isoprene. Our 
model is able to reproduce their vertical gradients, although the vertical profiles of 
HCHO and glyoxal appear to be sensitive to the choice of photochemical mechanisms. 
AM3BST underestimates HCHO by 32% and overestimates glyoxal almost by a factor of 
two near the surface. With the suppression of 5-ISOPO, in AM3B, HCHO is increased by 
17% and glyoxal is decreased by 38% near the surface, leading to better agreement with 


observations for both species. 


The abundance of glyoxal is also dependent on its heterogeneous loss. By assuming 
an irreversible reactive uptake on aerosols, an optimal value of Yglyx is selected to 
minimize the difference between modeled and observed glyoxal in the boundary layer. 
We find that the best agreement is achieved for AM3ST with yey, of 2 x 10° and for 


AM3B with yYelyx of zero (Table $2); we adopt those values in the following analysis. 
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Yelyx derived with AM3B is lower than the value from previous laboratory 
experiments (2.9 x 10° by Liggio et al. [2005b]), likely due to two reasons. First, there 
could be a missing source of gas-phase glyoxal in the current mechanism, which can be 
compensated by aerosol uptake and lead to a higher contribution to SOA formation. We 
find that approximate 30% increase of glyoxal production in AM3B would be required to 
compensate aerosol uptake of glyoxal with Ygiyx = 2.9 x 10°. Second, the effective aerosol 
uptake coefficient, as we assumed here to represent the net heterogeneous loss of glyoxal, 
is indeed smaller than the value from the previous study. One possibility is that the 
aerosol uptake of glyoxal is to some extent reversible [Galloway et al., 2009; Kroll et al., 
2005], which would lead to a lower estimate of yYgiyx. Another possibility is the 
suppression of glyoxal loss to ambient aerosols that are internally mixed with organic and 
inorganic components, resulting from its organic coating [Galloway et al., 2011b]. It is 
also possible that additional production of glyoxal in the particle-phase from other 
organic compounds reduces the net loss of glyoxal to the aerosol surface. In fact, the 
negligibly small uptake coefficients of glyoxal (0-8 x 10“ during the daytime and (2 + 1) 
x 107 during the nighttime) have also been determined from analysis of glyoxal 
observations in a previous field study in Los Angeles during the CalNex field campaign 
[Washenfelder et al., 2011}. 

Overall, the resulting average contribution of glyoxal to organic aerosols is about 0.8 
ug m* for AM3ST (Figure S7) and negligible for AM3B in the boundary layer over the 
Southeast U.S. We emphasize that this is likely the lower limit of glyoxal SOA, as other 
glyoxal sources such as anthropogenic VOCs and monoterpenes, are not included in our 


model. Accounting for these additional sources would require a higher sink of glyoxal, 
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and therefore results in higher glyoxal SOA. The estimate from AM3ST is comparable to 
previous model studies over the same region assuming reversible or irreversible reactive 
uptake of glyoxal onto aerosols [Knote et al., 2014; Li et al., 2015; Ying et al., 2015], 
while AM3B is not. We also find that the inclusion of IEPOX aerosol uptake may reduce 
glyoxal in the boundary layer by less than 15% over the Southeast U.S. (Figure S8). This 


would lead to an even lower estimate of glyoxal SOA. 
4.2 Rer and its dependence on NOx 


Comparison of modeled and observed Rep provides additional constraint on the gas- 
phase chemistry for glyoxal production. Figure 3 shows the co-variation of glyoxal and 
formaldehyde in the boundary layer during SENEX. The linear regression slope (Ror) 
derived from observations is 2.4% (R*=0.61) in the boundary layer over the Southeast US, 
suggesting that the dominant source of both glyoxal and formaldehyde in this region is 
isoprene [Kaiser et al., 2015]. The AM3ST mechanism overestimates Rep by a factor of 
two, while better agreement is achieved by the AM3B mechanism, due to the suppressed 
production of glyoxal and enhanced production of HCHO. Simulations with both 
mechanisms significantly underestimate the range of HCHO concentrations observed. 
Without heterogeneous loss of glyoxal, the model with AM3ST significantly 
overestimates Ror but predicts a similar result as assuming no heterogeneous loss of 


glyoxal with AM3B. 


We further examine the NO, dependence of HCHO, glyoxal and Rer during SENEX 
(Figure 4). A prominent feature of Figure 4 is the high similarity between the observed 


NO,-dependence of HCHO and glyoxal. Both HCHO and glyoxal are observed to 
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increase with NO, concentrations and to start to level off at 1 ppbv of NOx. Remarkably, 
Ror shows little variation across NO, concentrations, consistent with Kaiser et al. [2015]. 
Both AM3ST and AM3B are able to reproduce the NO, dependence of HCHO, although 
they tend to underestimate HCHO by 1-2 ppbv across NO, regimes, consistent with 
Wolfe et al. [2016]. AM3B shows better agreement with observations than AM3ST at 
high NO, level, indicating 5-channel is suppressed under ambient conditions. Using the 
optimized Ygiyx values described above (2 x 10° for AMBST and zero for AM3B), our 
model can roughly capture the NO, dependence of glyoxal. In contrast to Figure 1, 
neither AM3ST nor AM3B shows high glyoxal concentrations under low NO, conditions, 
reflecting slower production of glyoxal under such conditions due to low OH radical 
concentrations and slow ISOPO, production rate in the ambient atmosphere [Wolfe et al., 
2016]. Although the NO, dependent HCHO and glyoxal are individually captured 
qualitatively, a large positive bias in Rgris apparent for AM3ST across NO, regimes and 
for AM3B at NO, < 0.5 ppbv. We attribute the bias in AM3ST to the overestimate of 
glyoxal under low NO, and underestimate of HCHO across NO, regimes. The bias in 


AM3B is mainly due to the underestimate of HCHO. 


To examine the NO, dependence of HCHO, glyoxal and Re from MCM v3.3.1, the 
AM3M mechanism that mimics MCM v3.3.1 is then tested in our AM3 global model. 
Comparison to the SENEX data (Figure 4) shows that AM3M can also well reproduce the 
NO, dependence of HCHO, similar to the performance of other AM3 mechanisms. 
However, AM3M largely underestimates glyoxal and Rgr across all NO, levels, 
suggesting additional sources of glyoxal needed in the current MCM mechanism. 


Inclusion of heterogeneous loss of glyoxal would further reduce glyoxal in AM3M and 
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worsen the comparison. One possible reason for the underestimate of glyoxal in AM3M 
is the lack of a radical propagating channel for the IEPOXOO + HO: reaction in MCM 


v3.3.1, as described above. 


Ror provides a useful tool to estimate global production of glyoxal from isoprene 
oxidation. In particular, the fact that Ror is insensitive to NO, allows us to directly 
compute its global production rate. Assuming 0.4 HCHO yield per isoprene C (Figure 1, 
also consistent with Palmer et al. [2003]) and a global constant Reg of 2.4% (according to 
the measurements during SENEX), we can derive the global production of 23 Tg glyoxal 
with 500 Tg isoprene emitted annually [Guenther et al., 2006]. This is similar to a 


previous estimate of 21 Tg yr! glyoxal from isoprene [Fu et al., 2008]. 
5 Budget of Glyoxal in North America 


Table 1 summarizes the monthly averaged glyoxal budget over North America 
(20°~55° N, 60°~130° W) in the boundary layer (0-1.5 km) during June-July of 2013 for 
AMB3ST and AM3B. The total chemical production of glyoxal varies with the chemical 
mechanisms, from 0.32 Tg/month with AM3B to 0.44 Tg/month with AM3ST. The 
major sink of glyoxal with both mechanisms is photolysis, followed by aerosol uptake 
and OH oxidation for AM3ST and OH oxidation and aerosol uptake for AM3B, although 
the contribution of each pathway varies with mechanisms and yglyx. For example, in 
AMBST, heterogeneous loss is an important sink, accounting for 26% of the total 
chemical loss of glyoxal; however, the contribution is negligible in the AM3B 


simulations. Deposition accounts for only 1-2% of total glyoxal loss in the model. We 
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show in this work that the estimates of sources and sinks of glyoxal are dependent on the 


choice of gas-phase chemistry. 


6 Conclusions and Discussion 


A 0-D photochemical box model and a 3-D chemistry-climate model applied to data 
from the Southeast U.S. during the SENEX field campaign provide the first model 


evaluation of in situ glyoxal aircraft observations. 


We find that the three mechanisms (AM3ST, AM3B, and MCM v3.3.1) show 
similarity in HCHO but large differences in glyoxal, leading to opposite NO, dependence 
of Rer. Under low NO, conditions, the AM3 mechanisms predict much higher glyoxal 
yields than MCM v3.3.1, largely due to the significant contribution from isomerization of 
ISOPO, and oxidation of IEPOX via the ISOPO2+HO 2 pathway. Although the three 
mechanisms agree better under high NOx conditions, they show different pathways 


contributing to glyoxal production. 


With the constraints from field measurements during SENEX, AM3ST with an 
effective reactive uptake coefficient ygyx of 2 x 10° and AM3B without heterogeneous 
loss of glyoxal can best reproduce the observed vertical profile. The latter shows better 
agreement with observed Re in the boundary layer. These two choices lead to less than 
0.8 pg m> or negligible of glyoxal SOA in the boundary over the Southeast U.S., 
accounting for 0-14% of the total SOA in this region (Figure $7). These are likely the 
lower limit of glyoxal SOA due to missing sources of glyoxal that are not included in our 


model. Over North America, glyoxal sinks are dominated by photolysis, followed by 
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aerosol uptake and OH oxidation for AM3ST and OH oxidation and aerosol uptake for 


AM3B. Dry and wet deposition are negligible for both AM3 versions (Table 1). 


In the boundary layer, observations suggest a very similar NO, dependence for 
HCHO and glyoxal, resulting in a nearly constant Rgr across NO, levels [Kaiser et al., 
2015]. Both AM3ST and AM3B can roughly capture the NO, dependence of HCHO and 
glyoxal, although AM3ST tends to overestimate Rg p at all NO, regimes, likely due to the 
overestimate of glyoxal at low NO, and underestimate of HCHO across NO, regimes. 
AM3B shows a positive bias in Rgr at NOx < 0.5 ppbv due to insufficient production of 
HCHO. The MCM v3.3.1 like mechanism (“AM3M”) shows a large underestimate of 
glyoxal across all NO, levels. One possible reason is the lack of a radical propagating 


channel for the IEPOXOO + HO> reaction in MCM v3.3.1. 


While glyoxal has been studied extensively over the past decade to understand the 
magnitude of its heterogeneous uptake, we show here that its gas-phase production is a 
large source of uncertainty that requires equal consideration. We find that its production 
from isoprene oxidation varies greatly among different chemical mechanisms. This in 
turn greatly impacts global estimates of glyoxal and in particular its contribution to SOA, 
especially in regions with low to moderate NO, levels (Figure 1). Under high NOx 
conditions, models differ significantly in the production of glyoxal from oxidation of 
glycolaldehyde, as well as the fate of 65-[SOPO2. Under low NO, conditions, there is very 
little laboratory evidence available on the production of glyoxal from IEPOX, HPALDs, 
di-HPCARPs or other intermediate products. Although heterogeneous loss of IEPOX is 
not included in current work, sensitivity tests show that inclusion of this uptake reduces 


glyoxal concentrations by 10-20% for the mean vertical profile over the Southeast U.S. 
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Thus, the fate of IEPOX represents another mechanistic uncertainty for predicting 
glyoxal. Also, there are missing sources that are not represented in our model [Volkamer 
et al., 2015] would increase modeled glyoxal and thus also inferred SOA. Future 
laboratory measurements are urgently needed and may have important implications for 


understanding the contribution of glyoxal to SOA in past and future atmospheres. 
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832 Table 1. Budget of glyoxal over North America (20°~55° N, 60°~130° W) below 1.5 km 
833 in June-July of 2013. Percentage is the contribution of each pathway to the total chemical 
834 loss of glyoxal. 


AM3ST AM3B 
Yalyx (<10™) 2.0 0 
Production (Tg/month) 0.44 0.32 
Chemical Loss (Tg/month) 0.41 0.29 
Photolysis 57% 717% 
OH 17% 23% 
Uptake by Aerosols 26% = 
Dry Deposition (Tg/month) 1.2x107 1.2x107 
Wet Deposition (Tg/month) 1.1x107 1.6x107 
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Figure 1. Cumulative yields of glyoxal and HCHO in major pathways from isoprene 
oxidation at different NO, levels. Glyoxal and HCHO from isoprene oxidized by O3 and 
by NO3 are not shown due to low production in all the mechanisms. 
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Figure 2. Mean vertical profiles of isoprene, formaldehyde, NO,, glyoxal, organic 
aerosol (OA) and surface area density of aerosol during SENEX. Grey shades are the 
standard deviation (o) of the averaged profiles of the measured tracers. Dashed green line 
in panel (d) is model estimate without heterogeneous loss of glyoxal by AM3ST. 
Observed organic aerosol mass and aerosol surface area density are from dry particles. 
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Figure 3. Correlation between glyoxal and formaldehyde below 1.5 km during SENEX. 
Black dots represent observations, while green and red open circles show AM3ST with 
Yelyx = 2.0 X 10° (left) and AM3B without heterogeneous loss of glyoxal (right). Solid 
lines are linear regression lines, with regression slopes calculated from least-squares fit. 
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Figure 4. HCHO (ppbv), glyoxal (ppbv), and Rg in each NO, bin below 1.5 km during 
SENEX. Dots, dashed boxes and whiskers are the mean, interquartile range and lowest 
and highest of the observations; green squares and red triangles are the mean of AM3ST 
with Yglyx=2.0 x 10° and AM3B without heterogeneous loss of glyoxal respectively; blue 
dashed lines are the mean of AM3M without heterogeneous loss of glyoxal. 


